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ABSTRACT 

The phosphoesterase (PE) domain of the bacterial 
DNA repair enzyme LigD possesses distinctive 
manganese-dependent 3-phosphomonoesterase 
and 3 -phosphodiesterase activities. PE exemplifies 
a new family of DNA end-healing enzymes found 
in all phylogenetic domains. Here, we determined 
the structure of the PE domain of Pseudomonas 
aeruginosa LigD (PaePE) using solution NMR 
methodology. PaePE has a disordered N-terminus 
and a well-folded core that differs in instructive 
ways from the crystal structure of a PaePE»Mn 2+ » 
sulfate complex, especially at the active site that is 
found to be conformational^ dynamic. Chemical 
shift perturbations in the presence of primer- 
template duplexes with 3-deoxynucleotide, 
3 -deoxynucleotide 3 -phosphate, or 3 ribonucleotide 
termini reveal the surface used by PaePE to bind sub- 
strate DNA and suggest a more efficient engagement 
in the presence of a ^-ribonucleotide. Spectral per- 
turbations measured in the presence of weakly cata- 
lytic (Cd 2+ ) and inhibitory (Zn 2+ ) metals provide 
evidence for significant conformational changes at 
and near the active site, compared to the relatively 
modest changes elicited by Mn 2+ . 

INTRODUCTION 

Many bacteria have two independent mechanisms to 
repair DNA double-strand breaks (DSBs) — homologous 
recombination and non-homologous end joining (NHEJ) 
(1,2). DNA ligase D (LigD) is the central catalyst of bac- 
terial NHEJ in vivo (3). LigD is composed of three 



autonomous enzymatic modules: an ATP-dependent 
ligase (LIG) domain that joins 3'-OH and 5'-phosphate 
ends (4,5); a polymerase (POL) domain that adds 
ribonucleotides or deoxyribonucleotides to DSB ends 
and gaps (6-11); and a 3'-phosphoesterase (PE) domain 
(12) that hydrolyzes 3'-monophosphate termini to 3'-OH 
ends that can either be extended by POL or sealed by LIG. 
The PE domain also has a unique 3'-phosphodiesterase 
activity that resects short 3' ribonucleotide tracts created 
at DSB ends by the POL domain, to leave a 3'-OH 
monoribonucleotide end that is the preferred substrate 
for sealing by bacterial NHEJ ligases (12-16). 
Mutational studies of the PE domain of Pseudomonas 
aeruginosa LigD (PaePE) have identified residues essential 
for activity (13). Whereas His42, His48, Asp50, Arg52, 
His84 and Tyr88 are important for both the monoesterase 
and diesterase activities, Argl4, Aspl5, Glu21 and Glu82 
are implicated specifically in the monoesterase activity. 

In order to understand the structural basis for catalysis 
by the PE domain, we pursued crystallographic and 
solution NMR studies of PaePE in parallel. The 1.9- A 
crystal structure of the active phosphodiesterase core 
of PaePE revealed a novel fold in which an 8-strand 
P-barrel with a hydrophobic interior supports a distinctive 
crescent-shaped hydrophilic active site on its outer surface, 
wherein six of the essential side chains coordinate manga- 
nese and a sulfate anion (17). Subsequent crystal struc- 
tures of two archaeal PE proteins at high resolution — 
1.1 A for Candidatus Korarchaeum cryptofilm PE (CkoPE) 
and 2.1 A for Methanosarcina barkeri PE (MbaPE) — 
higlighted conservation of the .PaePE P-barrel fold, active 
site constituents and metal-binding site (18). Thus, 
PaePE exemplifies a novel family of 3'-end healing 
enzymes found in diverse taxa of the bacterial and 
archaeal domains of life. Putative PE homologs were 
also detected in several eukaryal proteomes (17). 
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Missing from the PaePE crystal structure was an 
N-terminal peptide segment, deletion of which selectively 
diminished the 3'-phosphomonoesterase activity of PaePE 
without affecting the 3'-phosphodiesterase activity (16). 
The MbaPE protein has the homologous N-terminal 
peptide, but electron density corresponding to this 
segment was missing in the MbaPE crystal structure, sug- 
gesting that it is disordered (18). CkoPE naturally lacks an 
equivalent of the N-peptide, and the associated 3'-phos- 
phomonoesterase activity. It was notable that the crystal- 
lization of PaePE, CkoPE and MbaPE could only be 
achieved in the presence of Mn 2+ and an oxyanion (either 
sulfate or phosphate), both of which were present as a 
Mn 2+ »oxyanion complex in the respective PE active sites 
(18). The phosphate and sulfate anions in the crystal struc- 
tures were taken to be mimetics of the scissile 3'-phosphate 
of the PE nucleic acid substrate. 

Our aim in the present study was to obtain a solution 
structure of PaePE in order to determine its conformation 
in the absence of metals or substrate mimetics and, if 
possible, to visualize the N-terminal peptide, in the absence 
of any perturbations including external non-physiological 
additives and crystal-packing forces. Initial NMR studies 
were performed with PaePE-(l-187), which produced 
spectra that were suitable for obtaining backbone assign- 
ments, but were not of high enough quality to yield a 
sufficient number of sidechain 'H assignments to allow 
structure determination. Serial C-terminal truncations 
identified PaePE-(l-177) as a catalytically active 3'-mono- 
esterase/3'-diesterase that produced spectra of consider- 
ably higher quality and yielded sufficient number of 'H 
assignments (19) suitable for calculating an NMR struc- 
ture of this highly dynamic enzymatic module. Here, 
we present the solution structure of PaePE-(l-177) 
determined using standard NOE-based approaches. By 
investigating spectral changes in the presence of each of 
three different oligonucleotide primer-template duplexes 
and a series of bivalent metal ions, we provide insight to 
the DNA binding and metal ion recognition properties of 
PaePE. 



MATERIALS AND METHODS 

Protein expression and purification 

Pc/ePE-(l-177) (referred to hereafter as PaePE) was pro- 
duced in E. coli as an N-terminal His 10 -Smt3 fusion; iso- 
tope labeling and purification of the tag-free PaePE 
protein are described in detail elsewhere (19). PaePE 
samples for NMR experiments were prepared in a buffer 
containing 50 mM Bis-Tris (pH 6.5), 200 mM NaCl and 
5mM DTT (referred to as the NMR buffer from hereon 
forward) and typically contained ~300uM protein. The 
lack of chemical shift changes in the 15 N,'H TROSY 
spectra acquired with and without a 10-fold excess of 
EDTA argued against the retention or acquisition of 
enzyme-bound bivalent metal ions during the purification 
process. Thus, all subsequent spectra were acquired in the 
absence of EDTA. 



Resonance assignment 

Backbone 'H, 15 N, 13 Ca, 13 C and sidechain 13 CP assign- 
ments were obtained via TROSY-based HNCO, 
HN(CA)CO, HNCACB and HN(CO)CACB experiments 
(20) utilizing uniformly 2 H,C, N-labeled PaePE. 
Sidechain 13 C and 'H assignments were obtained via 
(H)C(CO)NH and H(CCO)NH experiments utilizing uni- 
formly 13 C, 15 N-labeled PaePE (21). Experimental details 
have been described at length elsewhere (19). 

Distance and angular constraints 

Backbone 13 Ca, L3 Cp\ 13 C and 15 N chemical shifts were 
converted into dihedral <\> and v|/ angular restraints using 
the TALOS+ software suite (22). NOE-derived distance 
constraints to the backbone amide 'H nuclei were obtained 
using a 15 N-edited NOESY-HSQC (utilizing a uniformly 
15 N-labeled sample) with a 150-ms mixing time at a field 
strength of 800 MHz. Distance restraints involving side- 
chain protons were obtained from 3-dimensional 
13 C-edited NOESY experiments in D 2 0 (150 ms mixing 
time, 800 MHz), separately optimized for aliphatic and 
aromatic regions of the spectrum. Uniformly 
15 N, 13 C-labeled samples were utilized in both cases. 
NOESY crosspeak intensities were obtained using 
NMRView (23) and converted into distance restraints 
using the ambiguous distance restraints (ADR) protocol 
as implemented in the ARIA (24) software suite. Manual 
NOE assignments were supplemented with those assigned 
automatically using the ARIA scheme as described previ- 
ously (24). 

Structure calculations 

To facilitate structure calculation, the ordered regions of 
PaePE were determined using an improved pulse sequence 
to measure the 15 N-{'H} steady-state NOE at 600 MHz 
(25). Spectra were acquired with several repetitions of a 
T-Ti-i (x = 11ms) symmetric 'H saturation scheme for a 
total of 3 s. An additional delay of 2 s was used to allow 
lock stabilization. The reference experiment included a 5 s 
pre-delay. Spectra were acquired using sweep-widths of 
14ppm (512 complex points) and 32ppm (128 complex 
points) for the H and 15 N dimensions respectively, 
using 64 transients per t\ point. Based on the 5 N-{'H} 
steady-state NOE values, structure calculations were ini- 
tially performed excluding the flexible N- and C-termini of 
PaePE and including only residues Thr32-Glyl63. After 
the ordered 'core' had been suitably refined, the 
N-terminal 31 residues were included in the structure cal- 
culations to obtain the final structural ensemble. Few res- 
onance assignments were available for the disordered 
C-terminal region (164-177) that is dispensable for cata- 
lytic activity. Consequently, this polypeptide fragment was 
not included in the structure calculations. 

Manually assigned NOEs were combined with dihedral 
and hydrogen-bonding constraints into a standard assign- 
ment/structure-calculation protocol using the ARIA 2.3/ 
CNS 1.1 (24) software suite. Pro65 and Pro72 were con- 
strained to be in the cis conformation using the cw-proline 
patch during the structure calculation. Additional 
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assignments that were obtained automatically during a 
run were inspected for correctness and incorporated into 
subsequent runs. Structure calculations were performed 
using a Cartesian dynamics simulated annealing protocol 
that consisted of: (i) high temperature dynamics at 2000 K 
(10000 steps); (ii) 4000 steps of refinement; (hi) a 
Cartesian dynamics cooling stage from 2000 to 1000 K 
(6000 steps); and (iv) a second Cartesian dynamics 
cooling stage from 1000 to 50 K (5000 steps). Force con- 
stants of 5, 25 and 200kCalmol" 1 A" 2 for the dihedral 
constraints and 10, 10 and 50kCalmol _1 A -2 for the dis- 
tance constraints (ambiguous, unambiguous and hydrogen 
bond) were used during the three temperature stages of the 
Cartesian dynamics protocol. Prochiral Hp and Hy from 
long-chain amino acids and methyl groups from Val and 
Leu, were treated using the floating chirality method (24). 
For the final refinement step, a subset of the lowest energy 
structures were solvated in a 7.5 -A box of explicit solvent 
and refined using the full Lenard-Jones non-bonded po- 
tential and electrostatic interactions from the OPLS 
force-field using a protocol described in detail previously 
(26). In the final run, 2048 structures were calculated and 
the resultant 128 lowest energy structures were used in the 
water-refinement protocol. Fifteen of the lowest energy 
water-refined structures with the least number of viola- 
tions of experimental constraints were chosen to represent 
the final structural ensemble of PaePE (note that only resi- 
dues 1-163 were included in the final round of structure 
calculations). The quality of the resulting NMR structures 
was evaluated with PROCHECK (27). Regions of second- 
ary structure were defined using Promotif 3.0 (28,29). 

Preparation of oligonucleotides 

HPLC purified DNA oligos D16F1R1 (5'-TTTTTTGGG 
ACCGGTC(C 2 ' F )(rC); wherein the penultimate nucleo- 
side is 2'-fiuorocytidine and the 3' terminal nucleoside 
is ribocytidine), D18 (5'-TTTTTTGGGACCGGTCCC) 
and D18p (5'-TTTTTTGGGACCGGTCCCp: where p is 
a 3'-monophosphate) were purchased from Integrated 
DNA Technologies Incorporated. Stock solutions con- 
taining 100 ul of aqueous D16F1R1 (9.2mM) or 50 ul of 
aqueous D 1 8 ( 1 1 .6 mM) or 50 ul of aqueous D 1 8p (6 mM) 
were annealed in an Eppendorf thermocycler at 95°C for 
3min, at 60°C for 2min, at 40°C for lOmin and finally at 
25°C for lOmin, to generate duplex primer-templates each 
containing a 6-nt single-stranded tail. The concentrations 
of the annealed primer-template duplexes were estimated 
from absorption measurements at wavelengths of 260 and 
280 nm (see Figure 4a for details of the annealed 
duplexes). 

NMR measurements of protein-oligonucleotide 
interactions 

Samples of uniformly 2 H, 15 N-labeled PaePE (100 uM) in 
NMR buffer were used for NMR titrations. For all the 
primer-templates, 15 N,'H TROSY spectra were acquired 
for PaePE:oligo ratios of 1:0.05, 1:0.12, 1:0.2, 1:0.3, 
1:0.5, 1:1, 1:2 and 1:3 for D18; 1:0.05, 1:0.1, 1:0.2, 1:0.3, 
1:0.5, 1:1, 1:2 and 1:3 for D18p. Data were acquired 
at 600 MHz at 25° C using sweep-widths of 15ppm 



(512 complex points) and 36.2ppm (128 complex points) 
in the direct and indirect dimensions, respectively, using a 
recycle delay of 1.5 s and 48 transients per t\ point. For the 
D16F1R1 primer-template, the following PaePE:oligo 
ratios were used 1:0.05, 1:0.1, 1:0.2, 1:0.3, 1:0.5, 1:1, 1:2 
and 1:3 and 15 N,'H TROSY spectra were acquired 
at 800 MHz at 25°C using sweep-widths of 14ppm 
(512 complex points) and 32ppm (128 complex points) 
in the direct and indirect dimensions respectively using a 
recycle delay of 1.5 s and 64 transients per t\ point. 

Scaled chemical shift differences were calculated using 
the following equation: 

A8(i) = /(«„, H -5i,H) 2 +0.11(VN-M 2 (!) 

where <5 0 . H and <5 0 . N are the 'H and 15 N chemical shifts for 
a given residue in apo .PaePE, respectively, while S i u and 
<5 0 . N are the corresponding shifts at the zth titration point. 

NMR measurements of divalent cation binding 

Samples of uniformly 15 N, 2 H labeled PaePE (100 uM) in 
the NMR buffer (excluding DTT; removal of DTT 
produced no spectral perturbations) were used for the 
metal ion titrations (Mn 2+ , Cd 2+ and Zn 2+ ). Metal ions 
were added as dilutions from 7.5 mM stock solutions of 
MnCl 2 , CdCl 2 or ZnCl 2 . A series of 15 N, 'H HSQC 
datasets were obtained at 700 MHz (25°C) with sweep- 
widths of 14ppm (512 complex points) and 32ppm (128 
complex points) in the direct and indirect dimensions re- 
spectively, using recycle delays of 1.5 s (5.5 s for Mn 2+ ; to 
allow a more accurate comparison between the peak 
intensities in the reference state and those in the 
presence of the paramagnetic metal ion resulting in 
increased R\ rates and enhanced recovery toward equilib- 
rium between scans) and 32 transients per t\ point. The 
following protein-to-metal ratios were used: Cd + — 1:0.25, 
1:0.5, 1:1 and 1:2; Zn 2+ — 1:0.5, 1:1; Mn 2+ — 1:0.25, 1:0.5 
and 1:1. In addition, the effects of phosphate were 
analyzed using a 10-fold excess of phosphate in buffer 
containing equimolar ratios of PaePE and each of the 
three divalent cations. An additional dataset was collected 
at 600 MHz (25° C) for Zn 2+ in the presence of DTT using 
the following ratios of PaePE and Zn 2+ — 1:0.2, 1:0.4, 
1:0.6, 1:0.8, 1:1, 1:1.5 and 1:2. These data were fully con- 
sistent with those collected in the absence of DTT. The 
reversibility of the spectral perturbations induced by metal 
ion binding in all cases was tested by addition of an excess 
of EDTA. 

Fluorescence measurements of protein-DNA interactions 

A 1 uM sample of PaePE in the NMR buffer was used for 
the fluorescence measurements using a 3-ml cuvette that 
allowed continuous mixing. All measurements were per- 
formed on a Spex Fluorolog-3 (Horiba Jobin Yvon) at 
25°C. Tryptophan emission was monitored following ex- 
citation at 280 nm. Emission spectra were recorded from 
300 to 450 nm in 1 nm increments using a 0.5-s integration 
time and averaged over four scans. For the D16F1R1 
primer-template, the following PaePE:oligo ratios were 
used— 1:0, 1:0.05, 1:0.1, 1:0.2, 1:0.3, 1:0.5, 1:0.75, 1:1, 



1:1.5, 1:2, 1:3, 1:5, 1:8, 1:12, 1:18, 1:25 and 1:40. 
Experiments were performed first in the absence of diva- 
lent metal ions and then in the presence of 100 uM Mn 2+ . 
Apparent K d values were obtained by fitting the change in 
fluorescence intensity at the emission maximum in the 
absence of oligos (331 nm) with increasing D16F1R1 con- 
centration to a quadratic binding isotherm [Equation (2)] 
using in-house software that utilized the ODRPACK 
libraries (30). 

K d +[D]+[P] 0 -J{K d +[D]+[P] 0 ) 2 -4[D][P] Q 

Io ~ I - Io ° m> 

(2) 

where / is the intensity at 331 nm at a dsDNA concentra- 
tion of [D], I 0 is the corresponding intensity in the absence 
of DNA, is the intensity difference at saturation, [P] 0 is 
the protein concentration (held constant at 1 uM). 
Additionally, fluorescence titrations of 1 uM PaePE with 
Mn 2+ were performed in the absence of oligonucleotide 
using the following PaePE:metal ion ratios — 1:0, 1:0.05, 
1:0.1, 1:0.2, 1:0.5, 1:1, 1:2, 1:5, 1:10, 1:20, 1:40 and 1:80. 



RESULTS AND DISCUSSION 

Solution structure of PaePE 

Inspection of the final NMR ensemble of 15 structures (see 
Table 1 for statistics) reveals a disordered N-terminus 
(1-31). This fact is also borne out by the extremely low 
15 N-{'H} NOE values (0.08 ± 0.21; see Figure 1) for this 
region. However, some degree of local ordering is seen 
encompassing residues Aspl5-Thrl9 ( 15 N-{'H} NOE: 
0.18 ± 0.01). This region contains Aspl5 and is flanked 
by Argl4 and Glu21; these three sidechains are important 
for the 3'-monoesterase activity of PaePE (16). The distri- 
bution of c|)/\|/ angles for the Aspl5-Thrl9 segment indi- 
cates an overall extended structure (Supplementary Figure 
SI). A hydrogen bond between the carbonyl oxygen of 
Thrl9 and the amide of Glu21 is found in all 15 structures 
of the final ensemble. 

The core region (Gly33-Glyl63) of PaePE forms 
a P-barrel consisting of eight antiparallel P-stands 
(Figure 2) as follows: Leu35-His42 ((31), His48-Leu55 
(P2), Thr58-Ala63 (P3), Arg76-His84 (P4), Vail 05- 
Proll5 (P5), Hisl27-Glyl34 (P6), Serl38-Argl45 (P7) 
and Trpl55-Lysl59 (P8). A single a-helix 
(Proll9-Lysl25) bridges the P5 and p6 strands. The loca- 
tions of the secondary structural elements are consistent 
with those predicted from chemical shift data (19) and 
those seen in the crystal structure (in complex with 
Mn 2+ and sulfate) (17). A ninth P-strand, which packs 
against the N-terminus of P5 and includes residues 
Gly93 and Ser94, is seen in the majority of the structures 
in the solution ensemble; this strand extends to residues 
Phe91 and Ile95 in several of the structures in the 
ensemble. The crystal structure also shows a strand 
between residues Gly93 and Ile95. However, the low 
15 N-{'H} NOE values for this region (0.38 ± 0.15 for 
residues 91-95; Figure 1) and the elevated local 
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Total 


-4334.8 ± 106.7 


Van der Waals 


-303.8 ± 23.9 


Distance restraints 


321.5 ± 17.3 


RMSD from average structure^ 




All residues (236-298) 




Backbone (N, Cot, C) (A) 


0.46 ± 0.07 


Heavy atoms (A) 


1.29 ± 0.19 


Ramachandran statistics* 1 




Most favored region (%) 


65.0 (82.1) 


Additionally allowed (%) 


28.0 (17.2) 


Generously allowed (%) 


5.5 (0.6) 


Disallowed (%) 


1.5° (0.1) 



"Hydrogen bond restraints were H N -0 distance of 1.8-2.3 A and an 
N-O distance of 2.8-3.3 A. 

b Total number of violated constraints for the final ensemble of 15 
water-refined structures including residues 1-163. 
C RMSD (from the mean) for residues in regions of definite secondary 
structure 35^12, 48-55, 58-63, 76-84, 105-115, 119-125, 127-134, 
138-144, 155-159. 

d Ramachandran statistics for residues 33-163 for the 15 structures in 
the final water refined ensemble. The numbers in the parenthesis 
indicate the statistics for residues in regions of definite secondary struc- 
ture: 35—42, 48-55, 58-63, 76-84, 105-115, 119-125, 127-134, 138-144, 
155-159. 

e A majority of the residues found in the disallowed regions of the 
Ramachandran plot belong to Asp71 that flanks the cis Pro72. The 
remainder of the disallowed residues preceded or followed those for 
which no mainchain 'H, 15 N assignments were available. 



displacement in the NMR ensemble (Supplementary 
Figure S2) suggest its disorder in solution. This disordered 
strand ((39, Gly93-Ser94) is colored purple in Figure 2b. In 
addition, there are three instances of 3[ 0 helices, 
comprising residues Tyr88-Asp90 (in 40% of the struc- 
tures of the final ensemble), Gln97-His99 (20% of struc- 
tures) and Thrl46-Leul48 (20% of the structures). Only 
the first of these 3i 0 helices is seen in the crystal structure 
(17), albeit between residues Leu86-Asp90. There is no 
electron density in the crystal structure corresponding to 
the second 3 io helix, indicating significant disorder of this 
segment, as confirmed by the elevated spin-spin relaxation 
rates (i?2 at 900 MHz; data not shown) suggesting us-ms 
timescale dynamics. 

Analysis of the final NMR ensemble reveals high local 
displacement (Supplementary Figure S2) in the P-hairpin 
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Figure 1. 15 N-{'H} Steady-state NOE for PaePE. Data was acquired 
at 600 MHz. NOE values for the P-strands and the ot-helix are shown in 
red and blue respectively. 



between pi and (32 (Asp43-Tyr49, 0.69 ± 0.18 A, 
including the N-terminus of (32) and in the loop between 
(34 and (35 (Phe91-Aspl04, 0.40 ± 0.13A) that includes a 
dynamic 3i 0 helix as discussed above. Some disorder^ 
also seen in the loop between (37 and (38 (0.30 ± 0.09 A). 
Regions of definite secondary structure are generally well 
defined (0.10 ± 0.08 A) except the extreme N-terminus of 
(32. The regions of high disorder, i.e. elevated mean local 
displacement values in the loop regions, also display low 
^N-f/H} NOE values (Figure 1) indicating significant 
subnanosecond timescale dynamics. 

Differences between the NMR and crystal structures 
of PaePE 

Whereas the overall fold in solution is very similar to the 
crystal structure (17), significant differences can be seen, 
with the largest deviations being near the active site. The 
crystal structure (colored red in Figure 3) adopts a more 
closed conformation with a concave active site, compared 
to a more open conformation seen in the NMR ensemble 
in solution (the member of the NMR ensemble closest to 
the mean is colored cyan in Figure 3). The open-closed 
conformational switch is indicated by the dashed line with 
arrowheads in Figure 3. The largest distortions occur at 
the face of the P-barrel comprising strands pi, (32 and P5. 
Significantly smaller deviations are seen on the opposite 
face comprising strands P6, P7 and P8. In the crystal struc- 
ture, the catalytic Mn 2+ ion is coordinated with octahedral 
geometry by the sidechains of His42 (pi), His48 (P2) and 
Asp50 (P2), the sulfate anion and two waters. This 
network might favor a relatively rigid closed structure, 
whereas a more open structure, on average, forms in the 
absence of the metal ion and crystal packing forces. The 
C-terminal tail (177-187) in the crystal structure, missing 
from our PaePE NMR construct, tucks in under the pi- 
(32 P5 sheet and thus might also stabilize the closed con- 
formation seen in crystallo. Note that the NMR ensemble 
(Figure 2a) includes both open as well as closed struc- 
tures, though the former constitute a larger fraction of 
the ensemble. The archaeal PE crystal structures 




Figure 2. (A) NMR ensemble of 15 water-refined structures of PaePE 
(shown in cross-eyed stereo view) showing only residues Ser31-Glyl63. 
The p-strands and the a-helix are shown in red and blue, respectively. 
The disordered segment that lies between Phe91 and Asp 104 is shown 
in yellow (this includes a disordered 9th p-strand between Gly93 and 
Ser94). The P3— (34 loop containing the two proline residues in cis con- 
formation (Pro65 and Pro72) is shown in green. A 3 10 helix (indicated 
in cyan) forms in ~40% of structures in the final NMR ensemble. 
(B) Ribbon diagram for a representative structure (closest to the 
mean) of the NMR ensemble. A ninth p-strand that forms in 
majority of the structures between Gly93 and Ser94 is colored purple. 



(Supplementary Figure S3a) illustrate varying degrees of 
closure around the active site. In the present study, main- 
chain 'H, 15 N assignments could not be made for PaePE 
residues His42, Ala44, Ser45, Leu47, His48, Tyr49 and 
Tyr88 (though mainchain 13 Ca, 'Ha assignments could 
be obtained for all of these residues except His48). 15 N- 
{'H} NOE data for residues in this region where 'H, 15 N 
resonance assignments are available, suggest overall 
rigidity on the ps-ns timescale. However, substantial line 
broadening (confirmed by elevated R 2 values, data col- 
lected at 900 MHz; data not shown) associated with 
assigned resonances in this region indicate significant con- 
formational dynamics on thejis-ms timescale. The con- 
formational flexibility seen at the active site in the NMR 
ensemble is not evident in the crystal structure of the 
/V/ePE»Mn 2+ »sulfate complex. It therefore appears that 
the catalytic metal ion and the phosphomimetic sulfate 
are necessary to induce order at the active site and orient 
the important sidechains into their catalytically competent 
conformations. The extensive dynamics is also the likely 
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Figure 3. Differences between the NMR and X-Ray structures of PaeYE. The X-Ray structure (red) was solved in the presence of a catalytic Mr + 
ion and a phosphomimetic sulfate anion. The key residues that lie in the co-ordination sphere of the cation and anion are shown in green and blue, 
respectively. The regions where the differences are the largest are indicated — the (31— P2— (35 face near the catalytic site (dashed arrow), the (33-|34 loop 
(blue dashed oval) that contains the two cis prolines (trans in the X-Ray structure) and the (37— 138 loop (red dashed oval). Only the core region is 
shown in both cases. 



reason that PaePE (or any other PE domain) could not be 
crystallized in the absence of metal ions. 

Notable differences between the crystal and NMR 
structures are also seen in the long loop between (33 and 
(54 (circled in blue in Figure 3, right panel). Two proline 
residues (Pro65 and Pro72) in this segment were predicted 
to be in the cis conformation (albeit in a trans configur- 
ation in the crystal structure), based on the differences 
between their 13 CP and 13 Cy chemical shifts (19). 
Whereas the cis conformation for Pro65 could be con- 
firmed by the presence of Hot 1 " '-Ha 1 ) NOE cross-peaks 
(Supplementary Figure S4), extensive spectral overlap pre- 
vented the unambiguous assignment of either Ha 1 " -Ha 1 ) 
(characteristic of cis proline) or Ha^'-HS 1 (characteristic 
of trans proline) crosspeaks corresponding to Pro72. 
Nevertheless, the two prolines residues were locked in cis 
conformation for the structure calculations. The presence 
of cis prolines in the (33- P4 loop does not lead to signifi- 
cant deviations between the crystal and NMR structures 
for the strands P3 and P4 themselves, except for some 
minor differences at the extreme C- and N-terminal ends 
of P3 and P4, respectively. The biological role, if any, of 
the cis prolines is not immediately obvious though this 
loop contains two to four prolines in bacterial and 
archaeal PE domains (Supplementary Figure S5). 

The loop between P7 and P8 also displays differences 
between the crystal structure and the solution structure 
(circled in red in Figure 3). In the NMR ensemble, this 
loop folds back at an angle of ~90° about residues Asnl47 
and Serl53. This region is also quite dynamic on the ps-ns 
timescale with moderately low 15 N-{'H} NOE values 
(0.44 ± 0.09) compared with the flanking regions, a fact 
that is reflected by the relatively high mean local displace- 
ment values in the NMR ensemble (Supplementary Figure 
S2). Different degrees of closure for this loop are also seen 
in the case of MbaPE and CfcoPE (Supplementary Figure 
S3a). 

Comparison of the observed 13 Ca chemical shifts and 
those calculated from the crystal structures of PaePE 



using the program SHIFTX2 (31) averaged over three 
residues (Supplementary Figure S3b) reveals that largest 
differences occur around the regions that also display the 
largest structural deviations between the crystal and NMR 
structures. Most notable among these include the P1-P2 
hairpin and the P4-P5 loop. The largest residue-specific 
deviations occur for Leu86 (calculated to be shifted 
3.4 ppm downfield from that observed in solution) and 
Asp87 (3.13 ppm downfield) suggesting a higher degree 
of helicity in crystallo compared to that in solution, in 
agreement with that seen comparing the NMR and 
crystal structures (Figure 3). 

The sidechains of the catalytic residues in the NMR 
ensemble sample orientations similar to those seen in the 
crystal structure (Supplementary Figure S6). It is likely 
that the catalytically optimal orientations are locked in 
by the presence of the metal ion as discussed earlier. 
However, Glu82 (implicated in the monoesterase 
activity) samples sidechain orientations that are quite dif- 
ferent from those seen in the crystal structure. The reasons 
for this are not immediately apparent. Orientations for 
His48, for which no resonance assignments are available, 
were not calculated. 

In order to further probe the differences between the 
NMR and crystal structures of PaePE, we tested several 
alignment media (32) in an attempt to obtain a set of 
residual dipolar couplings. However, in spite of our best 
efforts, all these attempts failed because PaePE was 
unstable under all conditions tested. 

NMR analysis of the interactions of PaePE with DNA 

Positively charged residues distributed over the surface of 
PaePE (Supplementary Figure S7) could potentially act as 
loci for DNA binding. In order to determine the 
DNA-binding surface of PaePE, we analyzed spectral per- 
turbations induced in 15 N, ! H TROSY spectra of uniform- 
ly 15 N, 2 H-labeled protein by increasing concentrations of 
the 18-mer primer-template duplex, D18 (Figure 4a). 
Many spectral perturbations were seen in the presence of 
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Figure 4. (A) The three primer-template constructs used to measure 
spectral perturbations in PaePE: D18p contains a 3'-phosphate; 
D16F1R1. contains a 2'-fluorocytidine (to closely mimic a ribose) in 
the penultimate position and a 3'-terminal ribocytidine. (B) Results of 
the NMR titration of PaePE with D18. I5 N,'H TROSY spectra 
acquired at 600 MHz in the absence (red) and the presence of 
equimolar amounts of D18 (blue). The inset shows the trajectory of 
the resonance corresponding to Leu 157 that lies on p8 during the 
course of the titration. 



D18 (Figure 4b) reflecting the alteration of the chemical 
environment of the main chain nuclei due the combined 
effects of direct binding and induced conformational 
changes. Many residues were broadened beyond the 
threshold of detection early in the titration course 
(Supplementary Table SI). This extensive line broadening 
indicates exchange on the intermediate timescale (33). The 
broadened resonances correspond to residues situated 
near the active site, on the p3 P4 loop, on strands P4, 
(35 and (37, and on the (37-(38 loop. The positively 
charged residues Hisl27, Lysl36 and Argl40 were also 
completely broadened out early during the titration. 
From this group of perturbed residues, those on (35 
(Asp 109 and Glylll flank the positively charged 
ArgllO), and the set of positively charged residues that 
show localized perturbations (Lysl36, Argl40) are quite 
distant from the catalytic site, all lying on the opposite 
face of the protein (Figure 5c). If one or more of these 
residues were to take part in a productive binding event, a 
substantial distortion would be required of the duplex 
bound to this patch to access the catalytic site. A conform- 
ational change in this region in response to DNA binding 
at the opposite face, while also possible, is difficult to 



rationalize without further evidence. It is therefore likely 
that the spectral perturbations seen in the presence of the 
D18 primer-template result from both productive and 
non-productive binding events. In addition to these exten- 
sive line-broadening effects, significant chemical shift 
changes (Supplementary Figure S8) were also seen for 
several residues in and around strands P4 and (38. 

Given the ability of PaePE to catalyze the hydrolysis of 
the 3'-P0 4 of an all-DNA primer-template, we tested 
whether the presence of a phosphate group at the 3'-end 
(primer-template D18p) resulted in any alteration in the 
spectral perturbations. The overall distribution of the per- 
turbations (Supplementary Figure S8 and Supplementary 
Table SI) was quite similar to those seen for D18. Thus, 
the presence of a 3'-P0 4 group did not seem to have a 
substantial effect on the nature of spectral perturbations, 
at least in the absence of the catalytic Mn 2+ ions. 

The phosphodiesterase activity of PaePE requires a 
2'-OH in the penultimate nucleobase (12). Thus, successive 
base removal occurs until a single ribonucleotide remains 
at the 3'-end. Given this scenario, we tested the influence 
of a primer-template containing a 3'-terminal ribonucleo- 
tide on the spectral perturbations in PaePE. A primer- 
template construct (D16F1R1, Figure 4a) that was also 
2'-fluoro modified at the penultimate base (to mimic 
closely an RNA base) was used in the NMR titrations. 
In the presence of D16F1R1, substantial changes in the 
spectral perturbations were seen compared to D18 and 
D18p. The perturbations for residues comprising (35 
and the isolated perturbations (seen in the cases of the 
D18 and D18p constructs) near the positively charged 
residues Hisl27, Lysl36 and Argl40 were no longer seen 
(highlighted in Supplementary Table SI). While Lysl51 on 
the catalytic face was no longer significantly perturbed, the 
perturbation in Glyl50 persisted, albeit at a much higher 
concentration of DNA (Supplementary Table SI). Line 
broadening effects were limited to residues near the active 
site (Tyr37, Lys41 and Vail 05), the (33-(34 loop (Val64, 
Leu70), (34 (Arg76, Leu77, Ala78, Val79, Val81, Glu82 
and His84) and on the (37-P8 loop (Thrl46, Asnl47 
and Glyl50). In addition, large chemical shift changes 
(>0.11ppm at an equimolar ratio with D16F1R1, 
Figure 5d and Supplementary Figure S8) were seen near 
the catalytic site (Glu40: 0. 12 ppm), the |33-(34 loop (Val74: 
0.15 ppm; Lys75: 0.25 ppm), P4 (Gln80: 0.11 ppm), P7 
(He 144: 0.16 ppm), the P7-P8 loop (Glnl52: 0.22ppm; 
Serl53: 0.12ppm and Glnl54: 0.33 ppm) and P8 
(Trpl55: 0.37ppm; Phel56: 0.13ppm; Leul57: 0.2ppm 
and Lysl59: 0.15 ppm). Based on these observations, it 
appears that the binding site for DNA that is productive 
for catalysis involves a surface composed of portions of 
strands P3, P4, P7, p8 and intervening loops (Figure 5b 
and d). This region shows a high degree of conservation 
in PE domains found across the bacterial and archaeal 
taxa (Figure 5e). In the absence of a preferred ribose 
base at the 3'-end, detectable oligonucleotide binding 
seems to occur at additional putatively non-productive 
sites at the opposite face of the protein. The inability to 
fit the chemical shift or intensity changes with increasing 
oligonucleotide concentration to simple quadratic binding 
also suggests a more complex interaction dynamics than 
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Figure 5. (A) Residues important for the catalytic activity of PaeVE. Gln40 and Glu82 (along with Argl4 and Glul5 on the disordered N-terminus, 
not shown) important for the 3'-phosphomonoesterase activity are colored magenta; Lys66 and Arg76 that influence the 3'-phosphodiesterase activity 
are colored blue; His42, His48, Asp50, His84 and Trp88 influence both activities, are colored cyan. (B) Schematic representation of PaeVE showing 
the secondary structural elements in the same orientation (looking down into the catalytic cavity) as in (a) and the left panels of (C-E). The location 
of the catalytic residues as in (a) is indicated by the dotted oval. The strands comprising the p3— P4— 137— P8 face that shows the largest perturbations in 
the presence of oligonucleotide is indicated by the red dots, (c) Influence of D18 on PaePE. Residues that are broadened out beyond detection during 
the titration course are colored blue and residues that display chemical shift changes [A5, Equation (1)] >0.11ppm at an equimolar proteimoligo- 
nucleotide ratio, are colored red. The circled patches represent perturbations involving residues near pockets of positive charges on the face opposite 
to that bearing the catalytic residues. These perturbations are likely due to non-productive binding events. Residues for which data were not analyzed 
due to missing assignments, spectral overlap or weak peaks in the reference state, are colored green. Residues that show no significant spectral 
perturbations in the presence of oligonucleotide are shown in gray. The effects of D18p are quite similar to those of D18. (d) While the overall 
spectral perturbations on the catalytic face, in the presence of D16F1R1, are similar as in D18 and D18p, the perturbations on the face opposite to 
that bearing the catalytic residues, are not seen, (e) Sequence conservation in bacterial and archaeal PE domains are depicted on the PaePE surface 
using a cyan (least conserved) to maroon (most conserved) gradient. The catalytic residues and a large part of the putative oligonucleotide binding 
surface are well conserved. 
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a single-site binding event (also suggested by the 
fluorescence-based assays described below). It is tempting 
to speculate that the extensive line broadening effects 
seen could in part be due to exchange between the pro- 
ductive and non-productive sites. The presence of the 
ribose enables a higher degree of selectivity in the occu- 
pation of the productive sites. 

No significant changes in the chemical shifts were seen 
for the N-terminal segment with any of the primer- 
templates used in this study. No spectral perturbations 
were seen in the Gly98-Glyl01 segment (missing from 
the crystal structure); it was conjectured that this segment 
might grip the primer-template and position it at the cata- 
lytic site (18). 

No perturbations were detected for the N-terminal 
tail in a complex consisting of PaePE, D16F1R1 and 
Mn 2+ in a 1:0.5:0.5 ratio in a buffer containing 2mM 
phosphate (to mimic a state immediately following the 
3'-monoesterase activity in the presence of the product 
inorganic phosphate). Sub-stoichiometric amounts of 
the paramagnetic Mn 2+ and D16F1R1 were used to 
minimize the line-broadening effects due to the former 
and allow the visualization of a maximal number of res- 
onances (Supplementary Table SI). Based on the lack of 
observed perturbations at the N-terminus, it is conceivable 
that an intact phosphorylated 3'-end is necessary to 
observe the conformational changes, if any, required for 
the monoesterase activity. It is also possible that the con- 
formational changes are quite fast or transient on the 
chemical shift timescale and thus do not give rise to any 
observable chemical shift changes or line broadening 
effects. 

Fluorescence studies of the interactions of PaePE 
with DNA 

In addition to the NMR titration experiments discussed 
above, we measured native tryptophan fluorescence for 
PaePE to probe its interactions with D16F1R1. The 
protein contains five tryptophan residues (Trp62, 
Trpl08, Trpll3, Trpl41 and Trpl55), all of which are 
buried in the structure, as confirmed by a fluorescence 
emission maximum around 331 nm (34) following excita- 
tion at 280 nm. The change in fluorescence emission seen 
in the presence of D16F1R1 (Figure 6) was fit to a simple 
one-site binding model [Equation (2)] yielding an apparent 
K d value of 5.6 ± 0.6 uM. However, the quality of the fit 
to the single-site binding model is quite poor (Figure 6; 
black symbols), as evidenced by the nature of the re- 
siduals, clearly suggesting that this simple model is not 
suitable for PaePE/DNA interactions. Note, it is 
unlikely that this change in fluorescence emission is due 
to direct quenching by DNA, but rather due to a conform- 
ational change in the protein. Surface exposure of the 
tryptophan residues to allow contact with DNA would 
result in a red shift of the fluorescence maximum to 
around 340 nm (34). Based on the NMR titrations, the 
most perturbed tryptophan residues are likely to be 
Trp62, Trpl55 and Trpl58. 

Better fits to a quadratic isotherm were obtained in the 
presence of an excess of Mn 2+ (see Figure 6, red symbols) 
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Figure 6. Change in PaePE tryptophan fluorescence (measured at 
331 nm) in the presence of D16F1R1 with (red) and without (black) 
100 uM Mn 2+ . Experimental data are shown as filled circles and fits to 
a quadratic isotherm [Equation (2)] by solid lines. Almost no change in 
tryptophan fluorescence is seen with an excess of Mn 2+ ions alone (blue 
circles). The inset depicts the residuals for the fits in the absence (black 
dashed line) or the presence (red dashed line) of Mn 2+ . Note the 
reduced deviation from a simple quadratic binding isotherm and the 
improvement in the quality of the fits seen in the presence of Mn 2+ . 



with greatly improved residuals and a concomitant de- 
crease in the x 2 value per degree of freedom from 2509 
in the absence of Mn 2+ to 995 in its presence, albeit with a 
similar apparent K d value of 3.0 ± 0.1 uM. This suggests 
that the presence of the catalytic Mn 2+ ions in addition to 
a ribose base at the 3'-end further reduces binding at the 
non-productive site and helps localize the oligonucleotide 
substrate more efficiently near the active site at the pro- 
ductive binding surface. No significant change in fluores- 
cence emission was seen with Mn 2+ ions in the absence of 
DNA (see Figure 6; blue symbols). As stated earlier, 
PaePE requires Mn 2+ for catalytic activity but our 
results indicate that Mn 2+ alone does not induce conform- 
ational changes that are detectable by tryptophan fluores- 
cence. A comparison of the NMR ensemble with the 
crystal structure of Pc/ePE»Mn 2+ »sulfate complex indi- 
cates that the orientations of the tryptophan residues are 
quite similar. Therefore it is not surprising that the metal 
induced conformational changes are not reflected in our 
fluorescence data. However, as will be demonstrated in the 
following section, the overall conformational changes 
induced by the catalytic Mn 2+ are far more modest com- 
pared to those induced by Cd 2+ that partially supports 
catalysis and Zn 2+ that is highly inhibitory (12). 

Interactions with metal ions 

As mentioned previously, PaePE depends on Mn 2+ for its 
phosphomonoesterase and phosphodiesterase activities. 
Cd 2+ weakly supports the phosphodiesterase activity but 
not the subsequent hydrolysis of the 3'-P0 4 (12). In order 
to test whether Cd 2+ also engaged the residues that con- 
stitute the Mn 2+ binding site revealed by the crystal struc- 
ture, we analyzed spectral perturbations in PaePE in the 
presence of increasing amounts of Cd 2+ (Supplementary 
Figure S9). Resonances corresponding to almost all 
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residues (for which data could be analyzed) that lie on 
sheets pi, P2, p3 and p5 together with the intervening 
loops (see Figure 7b and Supplementary Figure S10) 
were progressively attenuated during the course of the ti- 
tration and finally broadened beyond the threshold of de- 
tection at an equimolar ratio of protein to Cd 2+ ions. 
Minimal further changes were seen upon increasing the 
metal ion concentration, suggesting a single ion-binding 
site on PaePE. The Cd 2+ -induced spectral perturbations 
could be reversed by addition of EDTA. Based on the 
nature of the spectral perturbations seen (no significant 
chemical shift changes or the appearance of new reson- 
ances), it appears that the PaePE-Cd 2+ interaction 



occurs on the intermediate exchange regime on the 
chemical shift timescale. 

Next, we tested the ability of PaePE to bind Zn 2+ that 
has been shown to be inhibitory in the presence of Mn 2+ 
(12). As in the case of Cd 2+ , strong spectral perturbations 
(Supplementary Figure S9), that reached their maximal 
value around a 1:1 stoichiometry, reversible by EDTA, 
were noted on the pi P2 P3 P5 face of the protein 
(Figure 7c). Additionally, some smaller but significant per- 
turbations, not present in the case of Cd 2+ ions (indicated 
by the black oval in Figure 7), were seen on the a-helix. 
More importantly, unlike in the case of Cd 2+ , a new set of 
resonances appeared during the course of the titration 




Figure 7. (A) Schematic representation of PaePE showing the secondary structural elements in the same orientation as in the left panels of (B-D). 
The largest conformational changes in the presence of metal ions involve the face formed by the strands pi, p2, f)3 and p5 (indicated by the red dots). 
Some weak binding is also seen in the a-helical segment [indicated by an asterisk, also circled in panels (b-d)]. Spectral perturbations in the presence 
of (b) Cd 2+ , (c) Zr? + and (d) Mn 2+ ions are mapped onto the PaePE surface. Reduction in peak intensity (I/Io) in 15 N,'H HSQC spectra of PaePE in 
the presence of an equimolar ratio of metal ions is depicted using a red (maximum attenuation in the presence of metal with respect to the reference 
state) to blue (minimum attenuation) scale (also see Supplementary Figure S10). Residues for which perturbations could not be analyzed due to 
missing assignments, spectral overlap or weak peaks in the reference state are shown in gray. Note that the signal attenuation in the presence of 
Mn~ + is due to a combination of conformational changes and paramagnetic relaxation enhancement effects. 
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(reaching their maximum intensity around an equimolar 
protein-to-metal ratio) suggesting an interaction in the 
slow exchange regime on the chemical shift timescale 
(see Supplementary Figure Sll for a representative 
example). Given the extensive nature of the spectral per- 
turbations seen in the case of both Cd 2+ and Zn 2+ , it 
appears they are the likely result of a conformational 
change involving the p 1 P2 (33 P5 face of the protein. 
This is also the region where the largest differences were 
seen between the NMR and crystal structures (Figure 3). 

Finally, we analyzed the spectral changes for PaePE in 
the presence of Mn 2+ (Supplementary Figure S9). Modest 
perturbations (Supplementary Figure S10), a combination 
of the effects of paramagnetic relaxation enhancement and 
conformational changes, were seen. The location of the 
line-broadening effects (see Figure 7d), though smaller in 
magnitude, were consistent with those seen in the cases of 
Cd + and Zn 2+ , suggesting, as expected, a common binding 
site and similar conformational changes in the presence of 
Mn 2+ . As in the case of Zn 2+ , additional line-broadening 
effects were also seen (Figure 7d, indicated by the black 
oval) localized on the oc-helical segment. It is therefore 
likely that the presence of several charged residues on 
this helix facilitate a weaker, non-specific interaction with 
bivalent metal ions. This binding is observable with Zn 2+ 
that appears to have the highest overall affinity toward 
PaePE and with Mn 2+ where the spectral perturbations 
due to binding events, even weak or transient ones, are 
greatly exaggerated due to relaxation enhancement in- 
duced by the paramagnetic nature of Mn 2+ . No significant 
differences in spectral perturbations were seen at an equi- 
molar protein-to-metal ratio in the presence or absence of 
a 10-fold excess of phosphate for any of the metal ions 
except perhaps in the case of Mn 2+ where some additional 
attenuation involving Lys41 and Asp43 (that flank the 
catalytic His42) was noted at the active site 
(Supplementary Figure S10). It is plausible that a 3'-P0 4 
facilitates the proper orientation of the catalytic Mn 2+ at 
the active site. Overall, our data suggest that whereas 
PaePE relies on Mn 2+ for catalysis, Mn 2+ induces only 
modest conformational changes compared to Cd 2+ or 
Zn 2+ . 

A crystal structure of CkoPE bound to Zn 2+ and phos- 
phate has recently been determined (35). The authors 
found that the structure of CA:oPE»Zn 2+ »phosphate 
complex is essentially identical to that of the previously 
determined CA:oPE»Mn 2+ »sulfate complex. However, the 
main difference lies in the coordination geometry of the 
metal ions: Zn 2+ and Mn 2+ were found to have tetrahedral 
and octahedral coordination geometries, respectively. This 
is in contrast to the large-scale spectral perturbations seen 
in our NMR data for PaePE in the presence of Zn 2+ (and 
to a certain extent Cd 2+ ) but not for Mn 2+ when compared 
to the metal-free enzyme. These differences cannot be ex- 
plained simply by an alteration of the co-ordination 
sphere of the metal, rather points to a significant conform- 
ational change involving a large part of the catalytic face. 
While it is possible that the conformational changes 
induced by Zn 2+ are different for the PaePE and CA:oPE 
domains [Zn 2+ acts as inhibitor for both proteins (12,35)], 



it is more likely that the conformers that display these 
changes are unstable under crystal packing forces. 



CONCLUSIONS 

We have obtained the solution structure of the PE domain 
of Pseudomonas aeruginosa LigD using solution NMR 
techniques. PE in solution displays the 8-stranded 
(3-barrel fold observed in the crystal structures of PaePE 
and two archaeal PEs. However significant differences 
between the solution and crystal structures are seen near 
the catalytic site. This involves strands (31, (32 and (35 that 
form a dynamic structure that samples multiple conform- 
ations on a (is-ms timescale leading to an open conform- 
ation, on average, compared to a more closed 
conformation seen in the crystal structure in the 
presence of Mn 2+ and sulfate. Additional differences are 
seen in the proline-rich (33-P4 loop where the two flanking 
proline residues are found to be in the cis conformation 
while trans conformations are seen in the crystal structure. 

Measurement of the spectral perturbations in the pres- 
ence of three different primer-template duplexes reveal the 
location of the DNA binding surface of PaePE to be at a 
large positively charged patch on the face that hosts the 
catalytic residues. In addition what appears to be non- 
productive binding of DNA is also seen around positively 
charged residues in two of the all-DNA primer-template 
duplexes. However this mode of binding is greatly reduced 
in the presence of a terminal ribonucleotide base that is the 
preferred substrate for LigD, as evident from the NMR 
spectral perturbations, and further decreased in the pres- 
ence of Mn 2+ ions, as revealed by fluorescence meas- 
urements. It therefore appears that a ribose 3'-end and 
the catalytic metal ion help localize the oligonucleotide 
efficiently near the catalytic site. In contrast, no 
DNA-induced perturbations are seen for the disordered 
N-terminus that has been implicated in the mononoeste- 
rase activity of LigD, whether in the absence or the 
presence of substoichiometric amounts of Mn 2+ . Thus, 
the nature of the contribution of this segment to the 
3'-monoesterase activity of LigD remains a mystery. 

Whereas Mn 2+ is essential for the catalytic activity of 
PaePE, we have found evidence for significant conform- 
ational changes in the presence of both Cd 2+ and Zn 2+ 
ions, with spectral features in the latter confirming the 
existence of an alternate conformation(s) in the presence 
of metal ions, this is in contrast to the minimal conform- 
ational changes seen by comparing the crystal structures 
of the C/foPE.Mn 2+ «sulfate and GtoPE«Zn 2+ «phosphate 
complexes, beyond the co-ordination sphere of the metal 
cation. Our spectral perturbation data opens up an 
intriguing possibility — it is plausible that within the con- 
formational ensemble sampled by PaePE in solution, there 
exists a dynamic equilibrium between conformers that are 
catalytically competent and others that are not. Zn 2+ ions 
help stabilize states compromised in their catalytic ability 
and greatly reduce the rate of exchange with the catalyt- 
ically competent states (or equivalently, shifts the equilib- 
rium toward the catalytically incompetent states) to an 
extent that this process is far slower than the catalytic 
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timescale leading to inhibition even in the presence of 
Mn 2+ ions. This would imply that the new resonances 
that appear in the presence of Zn 2+ correspond to an 
ensemble that is dominated by conformers unsuitable to 
induce catalysis. This shift in equilibrium is not as drastic 
in the case of Cd 2+ ions where the exchange between the 
conformers appears to be in the intermediate regime on 
the chemical shift timescale and is still fast enough to 
support some level of catalytic activity (12). 
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